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ABSTRACT 


A numerical implementation of shape design sensitivity analysis of 
built-up structures is presented, using the versatility and convenience 
of an existing finite element structural analysis code and its data base 
management system. This report is a continuation of the Technical 
Report 86-2, Part I: Conventional Design Parameters. The finite element 
code used in the implementation presented is the Engineering Analysis 
Language (EAL), which is based on a hybrid method of analysis. It has 
been shown that shape design sensitivity computations can be carried out 
using the database management system of EAL, without writing a separate 
program and a separate database. 

The material derivative concept of continuum mechanics and an 
adjoint variable method of design sensitivity analysis are used to 
derive shape design sensitivity informations of structural 
performances. A domain method of shape design sensitivity analysis and 
a design component method are used. Displacement and stress functionals 
are considered as performance criteria. 
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CHAPTER I 
INTRODUCTION 

1.1 Purpose 

To date there exist a wide variety of finite element structural 
analysis programs that are used as reliable tools for structural 
analysis. They give the designer pertinent information such as 
stresses, strains, and displacements of the structure being modeled. 
However, if this information reveals that the structure does not meet 
specified constraint requirements, the designer must make intuitive 
estimates on how to improve the design. If the structure is complex, it 
becomes very difficult to decide what step must be taken to improve the 
design. There is however, substantial literature [1] on the theory of 
shape design sensitivity analysis, which predicts the effect that 
structural shape design changes have on the performance of a structure. 

The purpose of this work is to develop and implement structural 
shape design sensitivity analysis, using the material derivative concept 
of continuum mechanics and the adjoint variable method [1] that takes 
advantage of the versatility and convenience of an existing finite 
element structural analysis code and its database management system. A 
domain method of shape design sensitivity analysis [2] is used in which 
design sensitivity information is expressed as domain integrals, instead 
of boundary integrals, to best utilize the basic character of the finite 
element method that gives accurate information not on the boundary but 
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in the domain. The finite element code that will be used is the 
Engineering Analysis Language EAL [3]. 

Using the full capabilities of the EAL system, design sensitivities 
can be calculated within the program, without knowing the source code of 
the program. This has the advantage that the user deals with only one 
program, with only one data base and no interfaces between different 
programs [4,5,6]. 

1.2 Continuum Approach of Design Sensitivity Analysis 

A number of methods could be used to implement structural shape 
design sensitivity analysis, but the most powerful is the continuum 
approach [1], This method can be implemented outside an existing finite 
element code [4,5,6], using only postprocessing data. This approach is 
convenient, because shape design sensitivity analysis software does not 
have to be embedded in an existing finite element code. The continuum 
approach to shape design sensitivity analysis calculation can also be 
implemented using a powerful database management system such as the 
Engineering Analysis Language (EAL) . Using the database management 
system of EAL, only one database is necessary for computation of shape 
design sensitivity information. That is, it is not necessary to crease 
interfaces and other datafiles to compute sensitivity information. 
Information on element shape functions used in the finite element model 
is, however, necessary for design sensitivity calculation. 

The continuum approach to shape design sensitivity analysis can 
easily be extended to complex structural systems that have more than one 
structural component using a design component method [7]. The shape 
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design sensitivity vector is the derivative of a constraint functional 
with respect to shape design parameters. The magnitude of each 
component reflects how sensitive the constraint functional is to a 
change in the associated shape design parameter. If a component of the 
vector is negative, the corresponding shape design parameter should be 
decreased to increase the value of the constraint functional. If a 
component of the vector is positive, the corresponding shape design 
parameter should be increased to increase the value of the constraint 
functional. In addition, if the magnitude of a component of the vector 
is large, then the corresponding shape design parameter will have a more 
substantial effect on design improvement. 

When a designer uses a finite element structural analysis code in 
design of a structure, it is most likely that a number of program runs 
are necessary before a substantially improved design is obtained. With 
the aide of a shape design sensitivity vector, the designer will know 
what direction to take to improve the design most efficiently. 



4 


CHAPTER II 

DESIGN SENSITIVITY ANALYSIS METHOD 

A detailed treatment of methods and procedures for calculating 
shape design sensitivities is given in Ref. 1, for constraint 
functionals such as compliance, displacement, stress and eigenvalues. 
For compliance and eigenvalue functionals, the adjoint equation is the 
same as the state equation, thus no adjoint equation needs to be 
solved. Each displacement and stress functional requires an adjoint 
load computation and an adjoint equation must be solved. Due to the 
symmetry of the energy bilinear forms, the state equation and the 
adjoint equation differ only in their load terms [1]. Using the reload 
or multi-load option of an existing finite element code, the adjoint 
equation can be solved efficiently [4] • For the displacement 
functional, the adjoint load is a unit load acting at the point where 
the displacement constraint is imposed. For the stress functional, the 
equivalent nodal force of the adjoint load has to be computed [1]* 

The flow chart of Fig. 1 shows the overall procedure. First, the 
model is defined by identifying the shape design parameters, constraint 
functionals, finite element model, and loadings. In the next step, EAL 
is used to obtain structural response. With the structural response 
obtained, an adjoint load is calculated for each constraint functional, 
external to EAL, using assumed displacement shape functions. The 
adjoint load is input to EAL, to obtain an adjoint response for each 
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Figure 1 . Flow Chart of Design Sensitivity Computation 
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constraint functional using reload option. Using the original 
structural response and the adjoint response, shape design sensitivity 
information is computed for each constraint, by numerically integrating 
the design sensitivity expressions. The process is convenient, since it 
uses data that are available in the database or easily computable 
outside EAL. 

The design sensitivity expressions contain derivatives of the 
displacement field. Because EAL uses a hybrid formulation for membrane 
and plate elements, a stress field is assumed and the displacement field 
is unknown. So derivatives of the displacement field cannot be computed 
directly. To overcome this difficulty, an acceptable displacement shape 
function is selected for the displacement field. Selection of the shape 
function is based only on the degree of freedom (nodal displacement) of 
the EAL finite element code. With nodal displacements calculated from 
EAL and using the selected shape function, derivatives of the 
displacement field at the Gauss points are calculated. An argument that 
supports this method is that, with the same degrees of freedom, 
different methods of finite element approximation give accurate results, 
if both approximation methods are acceptable, as in the case in 
contemporary FEM codes. 

To give the basic idea of implementation of shape design 
sensitivity analysis and computation procedures, a simple prototype 
structural component is investigated. Once shape design sensitivity 
analysis of a structural component is completed, the design component 
method of Ref. 6 can be used for design sensitivity analysis of built-up 


structures 
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Consider a thin elastic solid shown in Fig. 2, with thickness h(x) 
of the membrane, where x = [x^, X 2 ] T . The design is the shape of the 
domain of the membrane. The energy bilinear and load linear forms are 
[ 1 ], 


a (z,z) = // h(x) £ cr^(z) e^(z)dn 
u n i,j=l 


( 2 . 1 ) 


and 


£ (z) = J 7 T^z 1 dF 
U i=l 


( 2 . 2 ) 


where z = [z*, z^]^ is the displacement field, T = [T*, T^]^ is the 

i i i i — 

boundary traction, and a (z) and e J (z) are the stress and strain 
fields associated with the displacement z and the virtual displacement 
z, respectively. The variational state equation is [1] 


a (z,z) = £ (z) (2.3) 

u u 

for all kinematically admissible virtual displacements z. 

First consider the functional that represents the displacement z at 

A 

a discrete point x, 


= z(x) = // 6(x - x)z(x) dQ 

1 n 


(2.4) 


where 6(x) is the Dirac delta. The first variation of Eq. 2.4 is [1] 


= // S(x - x)z(x)dQ (2.5) 

ft 

The adjoint equation in this case is [1,3] 
a u (A,X) = // <5(x - x)X(x)dft 


( 2 . 6 ) 
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Figure 2. Clamped Plane Elastic Solid of Variable Thickness 

for all kinematically admissible displacements X. This equation has a 
unique solution X, which is the displacement field due to a unit point 

a 

load acting at a point x. Using the adjoint variable method, design 
sensitivity of the displacement functional is [1] 

♦I - JJo l [a ij (z)(7X i V ) + o^CX^Vz 1 V )]d0 
1 ” i,j-l J 3 

- /Jo [ I cr ij (z)e ij (X)]div Vdfl (2.7) 

i,j=l 


where X is the solution of Eq, 2.6, V is the gradient operator defined 
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as 7 = [ ] , V is the design velocity field and div V is the 

divergence of the design velocity field. To numerically integrate Eq. 
2.7, a two-by-two Gauss-point integration procedure is used. Using 
stress computation of membrane element E41 in EAL, stresses and strains 
can be expressed in matrix form as [3, 8] 

(o U (z), <r 22 (z), o 12 (z)) T . [p]{fi} (2.8) 

{c U (z), z 22 (z), e 12 (z)) T - [Er'M (2.9) 


where 




( 2 . 10 ) 


( 2 . 11 ) 


and [P] is the position coordinate matrix, which determines points where 
the stresses are obtained, {8} is the stress coefficient vector, and [E] 
is the elasticity matrix for a plane stress problem [9, 10]. To compute 
the derivatives of the displacement field z and the adjoint displacement 
field X which are required to evaluate Eq. 2.7, a bilinear displacement 
shape function is used. Because the membrane element E41 in EAL is a 
four noded element, a linear isoparametric element is assumed. 
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1 

u 


1 


V 



i 1 12 

vhere u and v , i = 1,2, 3, 4 are the nodal displacements of z and z , 

respectively. The individual shape functions are 


Nj = j (l - 5X1 - n) 

N 2 " J (1 + C)(1 - 
n 3 =|(1H)(1 + h) 

n 4 = J d - 5)(l + TO 


(2.13) 


The mapping of the finite element from the Xj-X 2 coordinate to 
the 5-h coordinate is given in Fig. 3. Using the chain rule of 
differentiation, the derivatives of the displacement field are 



N 

N 


1.5 

N 2 >; 

N 3,C 

N 4,5 

l.h 

N 2 >n 

N 3,n 

N 4 ,n 




(2.14) 


where [J] is the Jacobian of the mapping function from the 
coordinate to the 5 -T l coordinate. 
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Consider the general functional that represents a locally averaged 
stress on an element as 


% = // g(cr(z))m dfl 
Z fi * 


(2.15) 


where m p is a characteristic function, defined on a finite element ^ as 
— 


m P = { 


// dfl 
“p 


(2.16) 




V. 


and g is the stress function. The first variation of Eq. 2.15 is [1] 


% = III l ° lj ( z *) }® p dft 

Z a 1,3=1 3 o 1J * 

-If f [I g „(*> C ijkA (Vz kT V )J m dft (2.17) 

n i, j=l k,«rl o J 

+ // g(z)div V m dfi - // g(z)m dft // m div V dfi 

n p n p n p 

For the adjoint variable method, the material derivative z is replaced 
by a virtual displacement X in the first term on the right side of Eq. 
2.17, to define a load functional, and equate to the energy bilinear 
form to obtain the adjoint equation 


2 

a (X,I) = //[ l -%(z) a ij (X)] m dB (2.18) 

u a i , j = 1 3 a 3 

for all kinematically admissible displacements T. Equation 2.18 has a 
unique solution for a displacement field X [1], Using the adjoint 
variable method, design sensitivity of the stress functional is 
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V 2 


/Jo l [<> ij (z)(VX iT V .) + a ij (X)(Vz iT V )]dfl 

l,j=l 3 3 


- // n [ I a 1 ^(z)E i ^( X) ]div V dfi 

i,j=l 

- //„ I [Is li (*)c 13 kl (?. kT v.)J 

i ,3=1 k,£=l cr J 


(2.19) 


|mpd^ 


+ J/o g(z)dlv Vm D dfi - g(z)m p dfi ff Q m p div VdO 


i-j 

here C J is the stress-strain relation defined as 


ij „ijkfc k£ 
a J = C e 


i,j,k,Z = 1,2. 


( 2 . 20 ) 


To numerically integrate Eqs. 2.18 and 2.19, a two-by-two Gauss point 
integration procedure is used. 

If von Mises’ stress criteria is selected in the constraint 
functional, then 


g = [(a 11 / + (a 22 ) 2 * 3(a 12 ) 2 - o U o 22 } /2 


( 2 . 21 ) 


and 


■%- - C2o U - „ 22 )/2g 


do- 


ll 


iL. r,a 22 - o 11 


a o' 


22 


(2a - a )/2g 


iL_ = *J2 


3 a 


12 


= 3a /g 


( 2 . 22 ) 


which can be written in vector form as 


r_3g_i r dg 3g dg |T 
bo 11 ’ do 22 ’ do 12 * 


(2.23) 
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The equivalent nodal force is computed, based on the modified 
Hellinger-Reissner principle. The stress can be written as 

M = [P] {6} (2.24) 

(see also Eq. 2.8). The stress coefficients {b} can be expressed in 
terms of the nodal displacement coefficients {q } as [8] 

{6} = [h _1 ] [t] { q } (2.25) 

where H is the flexibility matrix. Using Eqs. 2.24 and 2.25, stress can 
be expressed in terms of the nodal displacements as 

{a} = [P] [H _1 ] [ T ] {q } (2.26) 

and thus the nodal equivalent forces for the adjoint load are 

IL { — ^fr } T [P] [H _1 ] [T]m dfi (2.27) 

A general structure is a collection of structural components that 
are interconnected by kinematic constraints at their boundaries. 

Results stated are from Refs. 1 and 7. The energy bilinear and load 
linear forms of a general system, consisting of membranes can be written 
as 


n - p 

a (z,z) = l [a (z,z) ] (2.28) 

< 1=1 

n . 

i (z) = l [f (z) ] 
u » , u 

fc=l 


(2.29) 
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where [a u (z,z)] and [ Jt^( z) ] are given in Eqs. 2.1 and 2.2. The state 
equation is [1] 


a u (z,z) = * u (z) 


for all kinematically admissible virtual displacement z. Since the 
energy bilinear and load linear forms of the state equation are just the 
sum of energy bilinear and load linear forms of each structural 
component, the design sensitivity equation of the system is a simple 
additive process [1, 7]. The generalized shape design sensitivity of a 
built-up structure for a displacement functional is 


n - I // 1 [a ij (z)(VX iT V ) + o ij (Vz iT V ,)]dfl 

^1 n * i,j-l 3 3 

- T If [ T a lj (z)e 1;i (X)]div V dft 

*=1 ^ i,j-l 


(2.30) 


where n is the number of components in the built-up structure. The 
shape design sensitivity of a built-up structure for a stress functional 
is 


n 2 T T 

r = I // , l [o lj (z)(7X 1 V.) + a i:i (X)(Vz i V ) ]dfl 
6 trl fl* i,j=l 3 3 

~ I If A I a 1- ^ (z) e 1 ^ ( X) ]div Vdfi 
fl i,j=l 

~ If a I [ l g i .(z)C ljk Vz k V )]m dfl 

i,j=l k,*Fl a J F 



(2.31) 
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where n is the number of components in the built-up structure and the 
stress is averaged over a subdomain of component q. 
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CHAPTER III 
PROGRAMMING ASPECTS 


So far analytical results and numerical algorithms for shape design 
sensitivity analysis have been stated. An outline of the basic 
organization of the EAL database management system is given in Part I of 
this report [6]. 

The program can handle two types of constraint functionals; 
displacement and stress. So far only membrane elements can be used to 
model a structure and evaluate shape design sensitivities. However, the 
method presented can be and will be extended to include truss, beam, 
plate, and three-dimensional elastic solid. A flow chart of the program 
is given in Figure 4. To use the program, the user sets the system 
control parameter, gives information about the design variables, 
specifies the constraints, sets up the finite element model and 
describes the velocity field. The finite element model is described in 
the runstream dataset INIT MODL 0 0, the velocity field is described in 
the runstream dataset VELO FILD 0 0, all other information is given in 
the runstream dataset PARA SET 0 0. After that, the program 
automatically computes sensitivities for the given constraints and 
design variables. System control parameters are 

LCAS - Actual load case number 

DE41 - Number of membrane elements 


NMDV 


- Number of design variables 
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CDIS - Number of displacement constraint functionals 

CS41 - Number of stress constraint functionals for element type 

E41 

CTOT - Total number of constraints 

For every constraint group, a table is required to describe the 
location of the constraint functional. For displacement constraints, 
two entries are needed for each constraint. The first entry is the node 
number on which the displacement sensitivity is evaluated and the second 
entry is the direction of the displacement constraint. For the stress 
constraints, a table is needed that gives the element numbers on which 
the stress constraint functionals are evaluated. The tables are 

CDIS List 00 - Displacement constraints 

ST41 List 00 - Stress constraints 

The velocity field is given in the runstream dataset VELO FILD 0 0. For 
each design parameter, one velocity field has to be given. The velocity 
for each nodal point has to be specified explicitly. 

Note that EAL stores all element information in an intrinsic 
coordinate system that depends on the displacement. Because derivatives 
of the displacement depend on the coordinate system used, the program 
recomputes element displacement in the element coordinates. 

The results of the shape design sensitivity analysis program in EAL 
are all stored in EAL-library file LI 2. If a stress constraint design 
sensitivity is specified in the input control parameter, the data set 
DVAL E41 1 0 is created to store the stress constraint functional 
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Set System 
Control Parameter 




[Give Initial FEM-Model 





Load 

Compute Sensitivity: 
Loop over all Elements 



Figure 4. Program Organization 












20 


values for constraints that are listed in the input data set ST41 
LIST. The shape design sensitivity parameter given in EAL-library file 
LI 2 have the following basic form: 

DSVE E41 "DV" "CIND” 

where 

DSVE E41 Design sensitivity parameter for membranes 

DV Design variable number 

CIND Constraint indicator 

between 21,000 and 22,000 Displacement constraint in 
global Xj-direction 

between 22,000 and 23,000 Displacement constraint in 
global X 2 ~direction 

between 23,000 and 24,000 Displacement constraint in 
global x^-direction 

between 40,000 and 50,000 Stress constraints in 
membrane elements 

For displacement and stress constraints, the last three digits of the 
constraint indicator "CIND" give the node number of the displacement 
constraint and the element number of the stress constraint, 
respectively. 
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Chapter IV 
Numerical Examples 

In order to check whether the shape design sensitivity information 
obtained is accurate, a comparison is made with the finite difference 
Ai|>. An appropriate design perturbation Au must be selected, in order to 
obtain a meaningful finite difference of the performance functional. 

That is, if Au is too small, At|> = ij> (u + Au) - t(u) may be inaccurate, 
due to loss of significant digits in the difference. On the other hand, 
if Au is too large, Aiji will contain nonlinear terms and the comparison 
with will be meaningless. 

To demonstrate the capability of the program, two examples are 
given in this chapter. First, a cubic box is presented in Section 
4.1. The cubic box is an extremely simplified model of a wing-box 
structure. A study of this square box can, however, provide a basis for 
study of the wing-box. The second example treated in Section 4.2, is a 
simple interface problem, composed of two plane elastic membranes with 
different material properties. 


4.1 Cubic Box Problem 

Consider a cubic box shown in Fig. 5. The box consists of five 
plane elastic components; top, bottom, two sides and end. Subdomains 
(or components) are numbered in Fig. 6 for convenience. The shape 
design variable of the system is height h of the box. A C°-velocity 
field is required on plane elastic component [11]. It suffices to use 



Figure 5. Cubic Box 



*2 r 02 


Figure 6. Numbering of Subdomains and Interfaces 
(view in (-) x-direction) 
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piecewise linear velocity fields on each component, which are given in 
Table 1. Note that Sh is design change. 

Table 1. Velocity Fields on Each Component 



External loads are applied along the edges of the top surface 
^13* **14* **15 w * t * 1 constant magnitude of 4.77 lb/ in in the positive Xj- 
direction. The cubic box is discretized by 320, 4-noded membrane 
elements (E41) having 377 nodes. Young's modulus and Poisson's ratio 
are 1.0 * 10^ lb/ in 2 and 0.316, respectively. The thickness of each 
component is 0.1 in. 

Several nodal points are selected to check accuracy of shape design 
sensitivity of the displacement functional of Eq. 2.7. Design sensitivity 
predictions ijij and finite differences = ^(u + Au) - (u) with 

6h = O.Olh are given in Table 2. 
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Table 2. Box Problem - Displacement Sensitivity 


Node 

Dir. 

♦ ,(u) 

* 10- 4 

4>.(u + Au) 
* 10' 4 

A*. 

* 10 -5 

* 10 3 

Ratio 

% 

129 

x 3 

0.7615 

0.7725 

0.1102 

0.1070 

97.1 

153 

x 3 

0.9449 

0.9544 

0.0951 

0.0917 

96.5 

157 

x 3 

0.7793 

0.7899 

0.1063 

0.1034 

97.3 

257 

4 

x 3 

0.3295 

0.3342 

0.0467 

0.0448 

96.0 

257 

1.4425 

1.4649 

0.2236 

0.2189 

97.9 

261 

A 

x i 

x i 

x 3 

x 3 

0.1448 

0.1467 

0.0195 

0.0190 

97.1 

277 

-0.1694 

-0.1713 

-0.0186 

*0.0179 

96.1 

281 

-0.3701 

-0.3745 

-0.0441 

*0.0435 

98.6 

281 

1.6099 

1.6302 

0.2024 

0.1995 

98.6 

285 

1.4964 

1.5184 

0.2201 

0.2148 

97.6 

To 

check 

accuracy of 

shape design sensitivity of 

stress 

functional 


of Eq. 2.19, 17 finite elements are selected. Design sensitivity 
predictions an< ^ finite differences AiJ^ with 6h * 0.01 h are given in 
Table 3. Both performance functionals yield excellent sensitivity 
results. 

For the computation of the shape design sensitivity analysis, the 
EAL runstream idea is used which is convenient, but the cost of direct 
disk read/write access and database overhead disk read/write access can 
be very high. The CPU-time and the overhead cost for the sensitivity 
analysis of the box problem are given in Table 4. 
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Table 

3. Box Problem 

- Stress 

Sensitivity 


El. No. 

* 2 (u) 

i|> 2 (u + Au) 

a* 2 

ih' 

v 2 

Ratio X 

1 

85.97 

87.14 

1.1725 

1.1598 

98.9 

2 

47.39 

48.01 

0.6157 

0.6288 

102.1 

9 

129.77 

131.48 

1.7100 

1.6932 

99.0 

17 

100.57 

101.64 

1.0660 

1.0544 

98.9 

18 

53.15 

53.74 

0.5911 

0.5842 

98.8 

25 

175.71 

177.16 

1.4540 

1.4418 

99.2 

65 

64.88 

65.76 

0.8741 

0.8644 

98.9 

81 

70.59 

71.41 

0.8143 

0.8046 

98.8 

129 

53.36 

54.08 

0.7224 

0.7142 

98.9 

130 

39.59 

40.13 

0.5377 

0.5316 

98.9 

145 

57.92 

58.62 

0.7064 

0.6980 

98.8 

146 

44.21 

44.72 

0.5143 

0.5081 

98.8 

225 

20.12 

20.38 

0.2648 

0.2621 

99.0 

248 

28.60 

28.85 

0.2503 

0.2474 

98.8 

249 

55.09 

55.38 

0.2960 

0.2919 

98.6 

256 

32.82 

33.20 

0.3796 

0.3756 

98.9 

284 

35.00 

34.94 

-0.0532 

-0.0533 

100.2 


Table 4. Overhead Cost 



FEM 

ANALYSIS 

Preprocessing 

Evaluation of each 
Performance Functional 

CPU- time 

474.8 

3791.6 

2375.4 

Total number of disk 

write accesses 

468 

50817 

44187 

Total number of disk 

read accesses 

1131 

240575 

166455 

Number of database 
overhead disk write 
access 

19 

9842 

2889 

Number of database 
overhead disk read 
access 

13 

145445 

67072 
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The overhead cost is given for three computational steps 
the finite element analysis 

- the preprocessing step that computes the stresses and 
displacement derivatives at the Gauss points for the 
actual load and prepares all necessary data for the 
sensitivity analysis. 

- the evaluation of each performance functional that 
evaluates the adjoint load, computes strain and 
displacement derivatives at the Gauss points for the 
particular adjoint load, and evaluates the sensitivity 
expression. 


4.2 Interface Problem 

Consider an interface problem, shown in Fig. 7. The structure 
consists of two plane elastic components, modeled with 45 nodes and 32 
membrane elements (E41) (see Fig. 8). The shape variation is 
parameterized by two design parameter bj and b£ that determine the 
position of the interface boundary and the height of the structure (see 
Table 5). 


Table 5. Velocity Fields 


Design Parameter Component 


Velocity Field 


1 1 

2 

2 1 
2 


V 1 - 0.1 Xjfibj 

V 1 - -0.1(x - 20) 6b. 

1 1 1 

V - 0.0 

v 1 - 0.0 


v 2 = 0.0 
v 2 - 0.0 

v 2 - 0.1(x 2 - 10)6b 2 
V 2 - 0.1(x 2 - 10)6b 2 


External load P = 100 N is applied to the structure at the upper left 
corner (x^ = 20.0, x 2 = 10.0). Young’s modulus and Poisson's ratio of the 






2 



Figure 7. Interface Model 
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structure are E » 2,0 * 
first component and E = 
the second component of 
0 • 1 cm • 

Several nodal points are selected to check accuracy of shape design 
sensitivity of displacement functional of Eq. 2.7. Design sensitivity 
predictions and finite differences Ai^ = i^(u + Au) - t|^(u) with 1% 
design perturbation are given in Table 6. 

Table 6. Interface Problem - Displacement Sensitivity 


10 5 N/cm^ and v * 0.316, respectively, for the 
1.0 x io 7 N/cm^ and r = 0.316, respectively, for 
the structure. The thickness for each member is 


Node 

No. 

Direction Design 
Var. 

♦i 

CM *H 

44 -. 

* io - 3 

♦;_ 3 
* 10 3 

Ratio 

X 

9 

1 

1 

0.0489747 

0.0480451 

- 0.9296 

- 0.9438 

101.5 

18 

1 

1 

0.0232287 

0.0227568 

- 0.4719 

- 0.4791 

101.5 

27 

1 

1 

- 0.000187702 

- 0.0001882 

- 0.0005 

- 0.000535 

107.0 

36 

1 

1 

- 0.023455 

- 0.0229872 

0.4678 

0.4748 

101.5 

45 

1 

1 

- 0.0469181 

- 0.045982 

0.9361 

0.9502 

101.5 

9 

2 

1 

0.162013 

0.157690 

- 4.323 

- 4.4083 

102.0 

18 

2 

1 

0.159341 

0.155000 

- 4.341 

- 4.4269 

102.0 

27 

2 

1 

0.158151 

0.153809 

- 4.342 

- 4.4275 

102.0 

36 

2 

1 

0.157767 

0.153427 

- 4.340 

- 4.4266 

102.0 

45 

2 

1 

0.157713 

0.153372 

- 4.341 

- 4.4264 

102.0 

9 

1 

2 

0.0489747 

0.0492151 

0.2404 

0.2437 

101.4 

18 

1 

2 

0.0232287 

0.0233417 

0.1130 

0.1147 

101.5 

27 

1 

2 

- 0.000187702 

- 0.000188328 

- 0.0006 

- 0.000533 

88.8 

36 

1 

2 

- 0.023455 

- 0.0235724 

- 0.1174 

- 0.1189 

101.3 

45 

1 

2 

- 0.0469181 

- 0.0471523 

- 0.2342 

- 0.2373 

101.3 

9 

2 

2 

0.162013 

0.162603 

0.590 

0.6025 

102.1 

18 

2 

2 

0.159341 

0.159913 

0.572 

0.5839 

102.1 

27 

2 

2 

0.158151 

0.158722 

0.571 

0.5834 

102.2 

36 

2 

2 

0.157767 

0.158339 

0.572 

0.5843 

102.2 

45 

2 

2 

0.157713 

0.158285 

0.572 

0.5845 

102.2 
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To check accuracy of shape design sensitivity of the stress 
functional of Eq. 2.19, 32 finite elements are selected. Design 
sensitivity predictions ^ and finite differences with 1% design 
perturbation are given in Table 7. Both performance functionals yield 
excellent results, with the exception of few stress sensitivity 
predictions, where the sensitivity is relatively small to other 
sensitivities and therefore small numerical errors in the overall 
calculation lead to large differences in the finite difference and the 
shape design sensitivity prediction. 

For the computation of the shape design sensitivity analysis, the 
EAL runstream idea is used which is convenient, but the cost of direct 
disk read/write accesses and database overhead disk read/write access 
can be very high. The CPU-time and the overhead cost are given in Table 
8 for the following computational steps; 

- Reading the runstream into the database. 

The finite element analysis. 

- The preprocessing step that computes the stress and 
displacement derivatives at the Gauss points for the actual 
load and prepares all necessary data for the sensitivity 
analysis. 

The evaluation of each performance functional for two design 
variables that evaluates the adjoint load, computes strain and 
displacement derivatives at the Gauss points for the 
particular adjoint load, and evaluates the sensitivity 
expression. 

- As an additional information, the overhead cost for the 
evaluation of each performance functional for one design 
variable is also given. 
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Table 7. Interface Problem - Stress Sensitivity 


Element No. 

Design 

variable 


A 


♦i 

Ratio 

X 

X 

1 

824.747 

808.933 

-15.754 

-16.071 

98.0 

1 

2 

824.747 

824.075 

- 0.672 

- 0.570 

117.9 

2 

1 

742.343 

727.787 

-14.556 

-14.778 

98.5 

2 

2 

742.343 

740.690 

- 1.653 

- 1.620 

102.0 

3 

1 

632.872 

620.581 

-12.291 

-12.479 

98.5 

3 

2 

632.872 

630.073 

- 2.799 

- 2.770 

101.0 

4 

1 

500.639 

490.968 

- 9.670 

- 9.817 

98.5 

4 

2 

500.639 

497.113 

- 3.526 

- 3.514 

100.3 

5 

1 

288.532 

283.359 

- 5.213 

- 5.257 

99.1 

5 

2 

288.532 

288.948 

0.416 

0.505 

82.3 

6 

1 

370.417 

364.384 

- 6.033 

- 6.120 

98.5 

6 

2 

370.417 

370.315 

- 0.102 

- 0.037 

27.3 

7 

1 

330.002 

325.300 

- 4.702 

- 4.759 

98.5 

7 

2 

330.002 

329.159 

- 0.843 

- 0.816 

103.4 

8 

1 

304.815 

300.636 

- 4.179 

- 4.228 

98.5 

8 

2 

304.815 

303.576 

- 1.248 

- 1.218 

102.5 

9 

1 

288.385 

283.201 

- 5.175 

- 5.267 

98.2 

9 

2 

288.385 

288.798 

0.413 

0.501 

82.5 

10 

1 

370.176 

364.152 

- 6.024 

- 6.112 

98.6 

10 

2 

370.176 

370.075 

- 0.101 

- 0.039 

26.2 

11 

1 

330.275 

325.668 

- 4.607 

- 4.665 

98.7 

11 

2 

330.275 

329.495 

- 0.780 

- 0.754 

103.4 

12 

1 

307.525 

303.585 

- 3.940 

- 3.990 

98.7 

12 

2 

307.525 

306.530 

- 0.995 

- 0.980 

101.1 

13 

1 

824.798 

808.992 

-15.806 

-16.061 

98.4 

13 

2 

824.798 

824.131 

- 0.667 

- 0.571 

116.8 

14 

1 

742.558 

728.017 

-14.541 

-14.762 

98.5 

14 

2 

742.558 

740.917 

- 1.641 

- 1.611 

101.9 

15 

1 

633.305 

621.021 

-12.284 

-12.470 

98.5 

15 

2 

633.305 

630.518 

- 2.787 

- 2.762 

100.9 

16 

1 

500.888 

491.189 

- 9.699 

- 9.845 

98.5 

16 

2 

500.888 

497.341 

- 3.547 

- 3.538 

100.3 

17 

1 

419.434 

411.272 

- 8.162 

- 8.284 

98.5 

17 

2 

419.434 

415.596 

- 3.838 

- 3.815 

100.6 

18 

1 

314.409 

308.864 

- 5.545 

- 5.628 

98.5 

18 

2 

314.409 

311.911 

- 2.498 

- 2.494 

100.2 

19 

1 

221.980 

218.903 

- 3.077 

- 3.105 

99.1 

19 

2 

221.980 

221.086 

- 0.894 

- 0.963 

92.8 

20 

1 

607.664 

605.998 

- 1.666 

- 1.860 

89.5 

20 

2 

607.664 

612.102 

4.438 

4.255 

104.3 

21 

1 

286.907 

282.791 

- 4.116 

- 4.166 

98.8 

21 

2 

286.907 

285.815 

- 1.092 

- 1.081 

101.6 

22 

1 

220.797 

218.224 

- 2.573 

- 2.618 

98.3 

22 

2 

220.797 

220.302 

- 0.495 

- 0.506 

97.8 

23 

1 

283.586 

282.316 

- 1.270 

- 1.278 

99.4 

23 

2 

283.586 

285.151 

1.565 

1.551 

100.9 

24 

1 

373.198 

370.467 

- 2.731 

- 3.148 

86.7 

24 

2 

373.198 

374.175 

0.977 

0.573 

170.4 

25 

1 

317.198 

314.061 

- 3.137 

- 3.179 

98.6 

25 

2 

317.198 

317.426 

0.228 

0.226 

101.1 
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Table 7 — Continued 


26 

1 

301.809 

299.335 

- 2.474 

- 

2.500 

98.9 

26 

2 

301.809 

302.232 

0.423 


0.424 

99.8 

27 

1 

295.251 

291.885 

- 3.366 

- 

3.388 

99.4 

27 

2 

295.251 

294.793 

- 0.458 

- 

0.445 

103.0 

28 

1 

184.710 

182.358 

- 2.352 

- 

2.391 

98.4 

28 

2 

184.710 

184.169 

- 0.541 

- 

0.543 

99.6 

29 

1 

421.246 

413.049 

- 8.197 

- 

8.322 

98.5 

29 

2 

421.246 

417.392 

- 3.854 

- 

3.838 

100.4 

30 

1 

294.280 

288.590 

- 5.690 

- 

5.775 

98.5 

30 

2 

294.280 

291.453 

- 2.827 

- 

2.822 

100.2 

31 

1 

169.748 

166.300 

- 3.448 

- 

3.499 

98.5 

31 

2 

169.748 

167.956 

- 1.792 

- 

1.793 

100.0 

32 

1 

51.074 

50.048 

- 1.026 

- 

1.043 

98.3 

32 

2 

51.074 

50.541 

- 0.533 

- 

0.534 

99.7 
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Table 8. Overhead Cost 



Reading in 
Program 

FEM 

Analysis 

Preprocessing 

Evaluation 
Each Performance 
Two Design 
Parameter 

of 

Functional 
One Design 
Parameter 

CPU- time 

24.4 

11.0 

252.2 

317.0 

198.8 

Total number of disk 
write accesses 

158 

78 

6133 

7689 

4817 

Total number of disk 
read accesses 

2 

152 

16,216 

22,994 

13,895 

Number of database 
overhead disk write 
access 

4 

14 

966 

469 

382 

Number of database 
overhead disk read 
access 

0 

14 

4391 

4591 

2949 
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Chapter V 
Conclusions 

Results of this study show that it is possible to combine the shape 
design sensitivity algorithms of Ref. 1 with a database management 
system of EAL. For the sensitivity computation, it is necessary to 
compute the derivatives of the displacement field and the adjoint 
displacement field. In the EAL finite element analysis, the stress and 
strain calculations are based on a hybrid formulation. However, assumed 
displacement shape functions are used for evaluating the displacement 
derivatives. Nevertheless, results of the shape design sensitivity 
analysis are very accurate, which indicates that it is not necessary 
that the sensitivity analysis and the finite element analysis use the 
same shape function. 

A database management system with a finite element capability and 
the adjoint variable method of design sensitivity analysis, permit 
implementation of a shape design sensitivity analysis method that does 
not require differentiation of element stiffness and mass matrices. It 
is shown that a database management system can be used to implement 
shape design sensitivity analysis, so only one program with one database 
is necessary. 
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